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Abstract. We construct a H^*-dynamicaI system describing the dynamics of a class of anharmonic 
quantum oscillator lattice systems in the thermodynamic limit. Our approach is based on recently 
proved Lieb-Robinson bounds for such systems on finite lattices [19] . 

1. Introduction 

The dynamics of a finite quantum system, i.e., one with a finite number of degrees of freedom 
described by a Hilbert space TC, is given by the Schrodinger equation. The Hamiltonian H is a 
densely defined self-adjoint operator on TC, and for a vector in the domain of H the state at 
time t satisfies 

(1.1) idtm = HiPit) . 

For all initial conditions ip{0) £ 7i, the unique solution is given by 

^(i) = e"'*-^V(0), for ah t G M. 

Due to Stone's Theorem e~^^^ is a strongly continuous one-parameter group of unitary operators 
on 7Y, and the self-adjointness of H is the necessary and sufficient condition for the existence of a 
unique continuous solution for all times. 

An alternative description of this dynamics is the so-called Heisenberg picture in which the 
time evolution is defined on the algebra of observables instead of the Hilbert space of states. The 
corresponding Heisenberg equation is 

(1.2) dtA{t) = i[H,A{t)], 

where, for each t S M, A{t) G B(7i) is a bounded linear operator on TC. Its solutions are given by a 
one-parameter group of *-automorphisms, Tf, of B{7i): 

A{t)=niAm)- 

For the description of physical systems we expect the Hamiltonian, H, to have some additional 
properties. E.g., for finite systems such as atoms or molecules, stability of the system requires 
that H is bounded from below. In this case, the infimum of the spectrum is expected to be an 
eigenvalue and is called the ground state energy. When the model Hamiltonian, H, is describing 
bulk matter rather than finite systems, we expect some additional properties. E.g., the stability 
of matter requires that the ground state energy has a lower bound proportional to A^, where is 
the number of degree of freedom. Much progress on this stability property has been made in the 
last several decades [Ml HI]. We also expect that the dynamics of local observables of bulk matter, 
or large systems in general, depends only on the local environment. Mathematically this is best 
expressed by the existence of the dynamics in the thermodynamic limit, i.e., in infinite volume. 
This is the question we address in this paper. 
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There are two settings that ahow one to prove a rich set of important physical properties of 
quantum dynamical systems, including infinite ones: the C* dynamical systems and the W* dy- 
namical systems p]. In both cases, the algebra of observables can be thought of a norm-closed 
*-subalgebra A of some algebra of the form B{7i), but in the case of the Vl^* -dynamical systems we 
additionally require that the algebra is closed for the weak operator topology, which makes it a von 
Neumann algebra. For a C*-dynamical system the group of automorphisms Tt is assumed to be 
strongly continuous, i.e., for all A £ A, the map 1 1-^ TtiA) is continuous in t for the operator norm 
(C*— norm) on A. In a VF*-dynamical system the continuity is with respect to the weak topology. 

In the case of lattice systems with a finite-dimensional Hilbert space of states associated with 
each lattice sites, such as quantum spin-lattice systems and lattice fermions, it has been known for 
a long time that under rather general conditions the dynamics can be described by a C* dynamical 
system, including in the thermodynamic limit When the Hilbert space at each site is infinite- 
dimensonal and the finite-system Hamiltonians are unbounded, this is no longer possible and the 
weak continuity becomes a natural assumption. 

The class of systems we will primarily focus on here are lattices of quantum oscillators but the 
underlying lattice structure is not essential for our method. Systems defined on suitable graphs, such 
as the systems considered in [6l |7] can also be analyzed with the same methods. In a recent preprint 
[T], it was shown that convergence of the dynamics in the thermodynamic limit can be obtained 
for a modified topology. Here, we follow a somewhat different approach. The main difference is 
that we study the thermodynamic limit of anharmonic perturbations of an infinite harmonic lattice 
system described by an explicit VF*-dynamical system. The more traditional way is to first define 
the dynamics of anharmonic systems in finite volume (which can be done by standard means |21j). 
and then to study the limit in which the volume tends to infinity. This is what is done in [1], but 
it appears that controlling the continuity of the limiting dynamics is more straightforward in our 
approach. In fact, we are able to show that the resulting dynamics for the class of anharmonic 
lattices we study is indeed weakly continuous, and we obtain a T4^*-dynamical system for the infinite 
system. The iy*-dynamical setting is obtained by considering the GNS representation of a ground 
state or thermal equilibrium state of the harmonic system. The ground states and thermal states 
are quasi-free states in the sense of [22], or convex mixtures of quasi-free states. In the ground state 
case the GNS representations are the well-known Fock reprensentations. For the thermal states the 
GNS representations have been constructed by Araki and Woods [2] . 

Common to both approaches, ours and the one of [1], is the crucial role played by an estimate of 
the speed of propagation of perturbations in the system, commonly referred to as Lieb-Robinson 
bounds [8l [TTl [16l [T71 [18] . Briefiy, if A and B are two observables of a spatially extended system, 
localized in regions X and Y of our graph, respectively, and Tt denotes the time evolution of the 
system then, a Lieb-Robinson bound is an estimate of the form 

\\[Tt{A),B]\\ < Ce-'^('^(^'^)-''l*l) , 

where C, a, and v are positive constants and d{X, Y) denotes the distance between X and Y. Lieb- 
Robinson bounds for anharmonic lattice systems were recently proved in [19] , and this work builds 
on the results obtained there. Our results are mainly limited to short-range interactions that are 
either bounded or unbounded perturbations of the harmonic interaction (linear springs). 

To conclude the introduction, let us mention that the same questions, the existence of the 
dynamics for infinite oscillator lattices, can and has been asked for classical systems. Two classic 
papers are [10^ I15j. Many properties of this classical infinite volume harmonic dynamics have been 
studied in detail e.g. [231 [9] and some recent progress on locality estimates for anharmonic systems 
is reported in [5| [20]. 

The paper is organized as follows. We begin with a section discussing bounded interactions. 
In this case, the existence of the dynamics follows by mimicking the proof valid in the context 
of quantum spins systems. Section 3 describes the infinite volume harmonic dynamics on general 
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graphs. It is motivated by an explicit example on Z'^. Next, in Section 4, we discuss finite volume 
perturbations of the infinite volume harmonic dynamics and prove that such systems satisfy a Lieb- 
Robinson bound. In Section 5 we demonstrate that the existence of the dynamics and its continuity 
follow from the Lieb-Robinson estimates established in the previous section. 

2. Bounded Interactions 

The goal of this section is to prove the existence of the dynamics for oscillator systems with 
bounded interactions. Since oscillator systems with bounded interactions can be treated as a 
special case of more general models with bounded interactions, we will use a slightly more general 
setup in this section, which we now introduce. 

We will denote by T the underlying structure on which our models will be defined. Here T will 
be an arbitrary set of sites equipped with a metric d. For T with countably infinite cardinality, we 
will need to assume that there exists a non-increasing function F : [0, oo) (0, oo) for which: 

i) F is uniformly integrable over F, i.e., 

(2.1) ||F|| := sup VF(d(x,y)) < oo, 
and 

ii) F satisfies 

^ , ^ F{d{x,z)) F{d{z,y)) ^ 

(2.2) C := sup > — — < oo. 

^Jrf^ F{d{x,y)) 

Given such a set F and a function F, by the triangle inequality, for any a > the function 

F,(x) =e-'^^F(x), 

also satisfies i) and ii) above with ||-Fa|| ^ 1 1 -^11 and Ca < C. 

In typical examples, one has that T dZ'^ for some integer d > 1, and the metric is just given by 
d{x,y) = \x—y\ = Yl'j=i\^j~yj\- In this case, the function F can be chosen as F(|x|) = (l-|-|x|)~'^~'^ 
for any e > 0. 

To each x G F, we will associate a Hilbert space Tix- In many relevant systems, one considers 
Tlx = i^^(M, dqx), but this is not essential. With any finite subset A C F, the Hilbert space of states 
over A is given by 

= (g)?^., 

and the local algebra of observables over A is then defined to be 

AA = <^i3{nx), 

xeA 

where B{7ix) denotes the algebra of bounded linear operators on Tix- 

If Ai C A2 , then there is a natural way of identifying C Aa^ ■, and we may thereby define 
the algebra of quasi-local observables by the inductive limit 

At = [j Ak, 

Acr 

where the union is over all finite subsets A C F; see [3111] for a discussion of these issues in general. 

The result discussed in this section corresponds to bounded perturbations of local self-adjoint 
Hamiltonians. We fix a collection of on-site local operators H^'^'^ = {Hx^x&v where each Hx is a self- 
adjoint operator over Tix- In addition, we will consider a general class of bounded perturbations. 
These are defined in terms of an interaction <I>, which is a map from the set of subsets of F to 
with the property that for each finite set X C F, $(X) G Ax and $(X)* = $(X). As with the 



B. NACHTERGAELE, B. SCHLEIN. R. SIMS, S.STARR, AND V. ZAGREBNOV 



Lieb-Robinson bound proven in [19], we will need a growth condition on the set of interactions $ 
for which we can prove the existence of the dynamics in the thermodynamic limit. This condition 
is expressed in terms of the following norm. For any a > 0, denote by Ba{T) the set of interactions 
for which 

Now, for a fixed sequence of local Hamiltonians H^°'^ = {Hx}xer, as described above, an in- 
teraction $ G Ba{T), and a finite subset A C F, we will consider self-adjoint Hamiltonians of the 
form 

(2.4) ifA = + Hf = ^Hx+ HX), 

xgA xca 

acting on TC\ (with domain given by 'S>xeA ^i^x) where D{Hx) C TLx denotes the domain of Hx)- 
As these operators are self-adjoint, they generate a dynamics, or time evolution, {t/^}, which is the 
one parameter group of automorphisms defined by 

T^{A) = e'^^^Ae-'^"^ for any A £ Aa- 

Theorem 2.1. Under the conditions stated above, for all t £M, AG Ar, the norm limit 

(2.5) ^mTHA)=n{A) 

exists in the sense of non- decreasing exhaustive sequences of finite volumes A and defines a group 
of *— automorphisms Tt on the completion of Ar- The convergence is uniform for t in a compact 
set. 



Proof. Let A C T be a finite set. Consider the unitary propagator 

(2.6) UA{t, s) = e**^^ e-^(*-^)^^A 
and its associated interaction-picture evolution defined by 

(2.7) rt%AA) = Ua{0, t) A UA{t, 0) for all A e Ar ■ 
Clearly, U\{t, t) = 1 for all t G M, and it is also easy to check that 

i^UA{t,s) = Hf{t)l(A{t,s) and - i^l(A{t,s) =UA{t,s) H'f{s) 
at as 

with the time-dependent generator 

(2.8) Hf{t) = e*^A^*//f e-^^A^* = ^ e^^A^* <I>{Z) e"*^A^* . 

ZcA 

Fix T > and X C T finite. For any A G Ax, we will show that for any non-decreasing, 
exhausting sequence {A„} of F, the sequence {t^^^{A)} is Cauchy in norm, uniformly for t G 
[— T, r]. Moreover, the bounds establishing the Cauchy property depend on A only through X and 
\\A\\. Since 



an analogous statement then immediately follows for {t^^"(A)}, since they are all also localized in 
X and have the same norm as 

Take n < m with X C A„ C Am and calculate 

(2.9) r;^™ (A) - T^%{A) = A {u^jo, s)UaAs, t) AUA,At, s)UaAs, 0)} ds . 
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A short calculation shows that 
d 



ds 



(2.10) 

where 
(2.11) 

and 



Ua^ (0, s) Uk„ {s, t) AUa^ (t, s) Ua^ {s, 0) 

iUA^{Q,s) [[Htjs) - HtSs)) MASs.t)AUAAt,s)\ Wa„(s,0) 



iZ^A„(0,.)e^^-A. B{s)y^-,[A{t) 
^(t) = e-*'^A;^e-"'"^ 



B{s) = e-^^t^(i/it(.)-i/X!(.))e 



(2.12) 



^ e**-^A^\A„$(2')e - ^ $(Z) 



ZcA 



ZCAn 



Jim \Jin 



ZCAm: 

ZnA,„\A„^0 



Combining the results of (j2.9p - (|2.12p . and using unitarity, we find that 



(2.13) 



< 



Ttt{A{t)],B{s) 



ds 



and by the Lieb- Robinson bound proven in [19j . it is clear that 



(2.14) 



< 



< 



T^l\ (Ait)) ,e'''-^Am\An$(Z)e"'^'^Am\An 



ZnAm\A„^0 



2PI 
2IIAI 



,2|l<I>|UCa|i-s| 



i) E E 11^(^)11 EE ^'^(^(^'^)) 

J/eA„i\A„ ZCAm: XSX ZS^ 



< 



2P||||$| 



< pi*IUC.|.-^l-l 5; ||<I>(Z)|| J]F„(<i(x,.)) 



y,z&Z 



» (e«"^"l*-l - 1) J] F,(d(x,z))F„(d(.,y)) 

J/GAm\A„ a^GX zGAm 



< 2PIIIICDIU (e2|l*ll"^«l*-^l-lj J] J]F,(d(x,y)). 

s/GA„\A„ xGX 

With the estimate above and the properties of the function Fa, it is clear that 



(2.15) 



sup 

tG[-T,T] 



as n, m — > 00 



and the rate of convergence only depends on the norm \\A\\ and the set X where A is supported. 
This proves the claim. □ 

If all local Hamiltonians are bounded, {tj } is strongly continuous. If the are allowed to be 
densely defined unbounded self-adjoint operators, we only have weak continuity and the dynamics 
is more naturally defined on a von Neumann algebra. This can be done when we have a suffiently 
nice invariant state for the model with only the on-site Hamiltonians. E.g., suppose that for each 
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X £ T, we have a normalized eigenvector 
define 



of Hx- Then, for all A £ A\, for any finite A C F, 



(2.16) 



xsA xeA 



p can be regarded as a state of the infinite system defined on the norm completion of ^r- The 
GNS Hilbert space Tip of p can be constructed as the closure of ^x^r V' G -^r (H)xer ^x- 

Then 



(2.17) 



{rt{A)-n,{A))' 



< 



(n(yl)-rf")(A))v. 
(rf "^(^) - ri^iA)) ^ + II (r(f")(^) - n^)) ^ 



For sufficiently large A„, the lim^^^g of middle term vanishes by Stone's theorem. The two other 
terms are handled by 12.51 It is clear how to extend the continuity to tp G Tip. 

We will discuss this type of situation in more detail in the next three sections where we consider 
models that include quadratic (unbounded) interactions as well. 



3. The Harmonic Lattice 

As noted in the introduction, we will consider anharmonic perturbations of infinite harmonic 
lattices. In this section we discuss the properties of the harmonic systems that we need to assume 
in general in order to study the perturbations in the thermodynamic limit. We will also show in 
detail that a standard harmonic lattice model posesses all the required properties. 



3.1. The CCR algebra of observables. We begin by introducing the OCR algebra on which 
the harmonic dynamics will be defined. Following [M], one can define the CCR algebra over any 
real linear space D equipped with a non-degenerate, symplectic bilinear form a, i.e. o" : P x P — > M 
with the property that if cr{f, g) = for all f £ "D, then g = 0, and 

(3.1) <y{f,g) = -(r{g,f) for all f,geV. 

In typical examples, V will be a complex inner product space associated with F, e.g. V = ^^(F) or 
a subspace thereof such as V = ^^(F), or ^^(Fq), with Fq C F, and 

(3.2) a{f,g) = lm[{f,g)] . 

The Weyl operators over D are defined by associating non-zero elements W{f) to each f £ T> which 
satisfy 

(3.3) W{f)* = W{-f) for each f eV, 
and 

(3.4) W{f)W{g) = e-^'^f^s^/^WU + g) for all /, G P . 

It is well-known that there is a unique, up to *-isomorphism, C*-algebra generated by these Weyl 
operators with the property that W{Q) = 1, W{f) is unitary for all / G and ||VF(/) — ]1|| = 2 
for all / G P\{0}, see e.g. Theorem 5.2.8 [4]. This algebra, commonly known as the CCR algebra, 
or Weyl algebra, over P, we will denote by W = W(P). 



ON THE EXISTENCE OF THE DYNAMICS 



7 



3.2. Quasi-free dynamics. The anharmonic dynamics we study in this paper will be defined as 
perturbations of harmonic, technically quasi-free, dynamics. A quasi-free dynamics on W(P) is a 
one-parameter group of *-automorphisms Tf of the form 

(3.5) TtiWif)) = W{Ttf), feV 

where : P — > P is a group of real-linear, symplectic transformations, i.e., 

(3.6) a{Ttf,Ttg) = a{f,g). 

As ||H^(/) — M^(5)|| = 2 for all / 7^ (7 € V, one should not expect rj to be strongly continuous; 
only a weaker form of continuity is present. This means that Tt does not define a C* -dynamical 
system on W, and thus we look for a T^*-dynamical setting in which the weaker form of continuity 
is naturally expressed. 

In the present context, it suffices to regard a W* -dynamical system as a pair {M, at} where M is 
a von Neumann algebra and at is a weakly continuous, one parameter group of *-automorphisms of 
M. For the harmonic systems we are considering, a specific -dynamical system arises as follows. 
Let p he a state on W and denote by {Hp,Trp,Qp) the corresponding GNS representation. We will 
assume that p is both regular and Tf-invariant. Recall that p is regular if and only if 1 1-^ p{W{tf)) 
is continuous for all / G P, and rj-invariance means 

(3.7) p{n{A)) = p{A) for all AgW. 

For the von Neumann algebra M, take the weak-closure of 7rp(W) in C{7ip) and let at be the 
weakly continuous, one parameter group of *-automorphisms of M obtained by lifting to 
The latter step is possible since p is r^-invariant, see e.g. Corollary 2.3.17 [3]. 

3.3. Lieb-Robinson bounds for harmonic lattices. To prove the existence of the dynamics 
for anharmonic models, we use that the unperturbed harmonic system satisfies a Lieb-Robinson 
bound. Such an estimate depends directly on properties of a and Tj. In fact, it is easy to calculate 
that 

[Tt{W{f)),W{g)] = {W{Ttf)-W{g)W{Ttf)W{-g)]W{g) 

(3.8) = [l-e'''^^^f^^^]w{Ttf)W{g), 



using the Weyl relations (j3.4p . For the examples we consider below, one can prove that for every 
a > 0, there exists positive numbers Ca and Va for which 



holds for all t G M and all f,g €z l'^{'L'^). In general, we will assume that the harmonic dynamics 
satisfies an estimate of this type. Namely, we suppose that there exists a number > for which 
given < o < oo, there are numbers Ca and Va for which 

(3.10) |l - e-(^*/'S)| < Cae'"^\'\ ^ \f{x)\ \g{y)\Fa {d{x,y)) 

holds for all t € M and all f,g £ ^^(r). Here we describe the spatial decay in T through the 
functions Fa as introduced in Section [2l Since the Weyl operators are unitary, the norm estimate 

(3.11) \\[rt{W{f)),W{g)]\\<Cae"^\'\Y.\f{x)\\g{y)\Fa{d{x,y)) , 

x,y 

readily follows. 
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3.4. An important example. Using the example given below, we illustrate the general discussion 
above in terms of a standard harmonic model defined over T = U^. We begin with a description 
of some well known calculations that are valid for these models when restricted to a finite volume. 
This analysis motivates the definition of the harmonic dynamics in the infinite volume. We then 
demonstrate that this infinite volume dynamics satisfies a Lieb-Robinson bound. By representing 
this dynamics in a suitable state, the relevant weak-continuity is readily verified. Interestingly, our 
analysis also applies to the massless case of w = 0, see below, and we discuss this briefiy. We end 
this subsection with some final comments. 

3.4.1. Finite volume analysis. We consider a system of coupled harmonic oscillators restricted to a 
finite volume. Specifically on cubic subsets = (— L,L]^ C Z,'^, we analyze Hamiltonians of the 
form 

d 

(3.12) Hi= j2pI + ^\i + Y1 - 

acting in the Hilbert space 

(3.13) Ha^ = (g) L\R,dq,). 

Here the quantities px and Qx, which appear in (j3.12p above, are the single site momentum and 
position operators regarded as operators on the full Hilbert space TCa^ by setting 

(3.14) = 110 • • • (g) 1 (g) -i-^ (g) 1- • • (Xi 1 and g^ = ]l(g)---(gl(gg(gl---(g]l, 

dq 

i.e., these operators act non-trivially only in the x-th factor of Ha^- These operators satisfy the 
canonical commutation relations 

(3.15) \Px,Py] = [qx,qy] = and [qx,Py] = iSx,y, 

valid for all x,y ^ K^. In addition, {ej}^^^ are the canonical basis vectors in l/, the numbers 
\j > and w > are the parameters of the system, and the Hamiltonian is assumed to have 
periodic boundary conditions, in the sense that qx+ej = (lx-{2L~i)ej if x E but x + ej A^. It 
is well-known that Hamiltonians of this form can be diagonalized in Fourier space. We review this 
quickly to establish some notation and refer the interested reader to [19] for more details. 
Introducing the operators 

(3.16) Qk = = E 



defined for each k G A^ = : xGA^}, and setting 



W2 



+ 4^Aj sm'^{kj/2) 



(3.17) j{k) 
one finds that 

(3.18) = E ^(^^) i'^Kh + i) 



keAi 



where the operators bk and bt satisfy 



(3.19) bk = J—^ Pk - i\ ^ Qk and bl = P^k + Q~k 

In this sense, we regard the Hamiltonian as diagonalizable. 
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Using the above diagonalization, one can determine the action of the dynamics corresponding to 
on the Weyl algebra >V(^^(Al)). In fact, by setting 



(3.20) 



W{f) = exp 



Re[/(x)]g. + Im[/(x)]p. 



for each / G £^(Al), it is easy to verify that (13. 3|) and (13. 4p hold with a-{f,g) = lm[{f,g)]. It is 
convenient to express these Weyl operators in terms of annihilation and creation operators, i.e., 



(3.21) 



{Qx + iPx) and a* = {q^, - ip^) 



which satisfy 

(3.22) [ax, ay] = [a* , a*] = and [ax,a*] = 5x,y for ah x,y £ Al 

One finds that 

i 



(3.23) T^(/) = exp _ (a(/) + a*(/)) 
where, for each / G £'^{Al), we have set 

(3.24) a{f) = ^ 7M«x, a*if) = ^ fi 



xj a^ 



x£Ar 



Now, the dynamics corresponding to if^, which we denote by r/", is trivial with respect to the 
diagonalizing variables, i.e., 

(3.25) rHbk) = e-^'^^'^Hk and r^^ (bl) = e^'^^l , 
where 6^ and are as defined in (|3.19p . Hence, if we further introduce 

(3.26) 6, = 



e^'^^bk and 



bi = Yl ^^'''^^ 



for each x £ Al and, analogously to (j3.24p . define 

(3.27) bif) = Y7{x)bx, h*{f) = Y f{x)bl, 

xgAl x£Al 

for each / S £^(Ai), then one has that 

(3.28) Ti^ib{f)) = b{[T-HltJ']f) , 

where J- is the unitary Fourier transform on ^^(A^) and Mt is the operator of multiplication by 
g2i7(fc)f -j^ Fourier space with ^{k) as in (j3.17p . We need only determine the relation between the 
a's and the 6's. 

A short calculation shows that there exists a linear mapping U : £^(Ai) — > ^^(A^,) and an 
anti-linear mapping V : £'^{Al) — > £^(Ai) for which 

(3.29) bif)=aiUf) + a*iVf), 

a relation know in the literature as a Bogoliubov transformation [13]. In fact, one has that 
(3.30) 



U = -/"^Ufr^J^ and V = -T~^Mr TJ 
2 + 2 



where J is complex conjugation and M^j^ is the operator of multiplication by 
(3.31) r±(A;) = -^±y^, 
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with j{k) as in (j3.17p . Using the fact that T± is real valued and even, it is easy to check that 

(3.32) U*U -V*V = 1 = UU* -VV* 
and 

(3.33) V*U -U*V = = VU* - UV* 

where we stress that V* is the adjoint of the anti-linear mapping V. The relation (j3.29p is invertible, 
in fact, 

(3.34) 

and therefore 



aif)=biU*f)-b*iV*f), 



(3.35) 



W{f) = exp 



V2 



{b{{U* - V*)f) + b*{{U* - V*)f)) 



Tt{W{f)) = W{Ttf) , 



Clearly then, 
(3.36) 

where the mapping Tt is given by 

(3.37) Tt = {U + V)J'^^Mtr{U* - V*) , 
and we have used ()3.28p . 

3.4.2. Infinite volume dynamics. It is now clear how to define the infinite volume harmonic dynam- 
ics. Consider a subspace V C £^(Z'^) and define W{V) as above with (j{f,g) = Im[(/, 5)]. First 
assume u> > 0, take 7 : [— 7r,7r)'^ — > M as in (j3.17p . and set U and V as in (j3.30p with (j3.3ip . If 
u; > 0, both U and V are bounded transformations on £^(Z'^). We will treat the case a; = by a 
limiting argument. The mapping Tt defined by setting 

(3.38) Tt = {U + V)J^^^MtJ^{U* - V*) , 

is well-defined on £^(Z'^). To define the dynamics on yV('D) we will need to choose subspaces T> 
that are Tt invariant. On such P, Tt is clearly real-linear. With (j3.32p and (j3.33p . one can easily 
verify the group properties Tq = 1, T^j^t = Tg oTt, and 

(3.39) lm[{Ttf,Ttg)]=lm[{f,g)] , 

i.e. Tt is sympletic in the sense of (j3.6p . Using Theorem 5.2.8 of [3], there is a unique one parameter 
group of *-automorphisms on W(T>), which we will denote by Tt, that satisfies 

(3.40) Tt{W{f)) = W{Ttf) for all /G P. 

This defines the harmonic dynamics on W{T>). 

Here it is important that Tt : T> ^ T>. As we demonstrated in [19], the mapping Tt can be 
expressed convolution. In fact, 



(3.41) 
where 

(3.42) 



1 



1 

(27r)<^ 
1 



Im 
Re 
Im 



1 



(-1). 



7(fc) e*(^-^-2T(fc)*) dk 



Using analysis similar to what is proven in [19], the following result holds. 
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Lemma 3.1. Consider the functions defined in ^3.4^ - For uj > 0, Ai, . . . , > 0, but such that 
C(^,A = ("^^ + '^Yl'j=i '^jY^'^ > 0' '^^^ '^''^y ^ > 0; bounds 

^ ^-M(k|-c..Amax(|,e("/2)+l)|i|) 



(3.43) 



^ ^-/.(ixl-c^,;, max(| , eW2)+i) \t\) 



hold for all t and x G Z"'. Here \x\ = Yl'j=i 

Given the estimates in Lemma l3.ll equation (|3.4ip and Young's inequality imply that Tf can be 
defined as a transformation of £^(Z'^), for p > 1. However, the symplectic form limits us to consider 
V = iP{Z'^) with 1 < p < 2. 

The following bound now readily follows: 



\lm{Ttf,g)\ < (l + 26^^/2^^^^ + 2c-\^ 
(^•44) x^|/(.)|b(,)|e-Kl-l--.^-(f-'''^^'^Ol*l) 

This implies an estimate of the form p.9p . and hence a Lieb-Robinson bound as in (|3.1ip . 
A simple corollary of Lemma |3. II follows. 

Corollary 3.2. Consider the functions defined in \3.4S^ - For cu > 0, Ai,...,Ad > 0, but with 
Cu),\ = (t^^ + '^^'j=i ^jY^"^ > 0- Take || • ||i to 6e the i^-norm. One has that 

(3.45) -Jolli ^0 as t ^ 0, 
and 

(3.46) ||if^"^||i^O as t^O, form £ {-1,1}. 

Proof. The estimates in Lemma |3 . 1 1 implv that the functions H^"^^ are bounded by exponentially 
decaying functions (in \x\). These estimates are uniform for t in compact sets, e.g. t G [—1, 1], and 
therefore dominated convergence applies. It is clear that H^\x) = 6o{x) while Hq"^\x) = for 
m £ { — 1, 1}. This proves the corollary. □ 

3.4.3. Representing the dynamics. The infinite-volume ground state of the model ()3.12p is the 
vacuum state for the 6— operators, as can be seen from ()3.18p . This state is defined on W{T>) by 

(3.47) p(H^(/)) = e-3ll(^*-^*)^ll' 

By standard arguments this defines a state on W(V) [i]. Using (j3.38p . (|3.32p and (|3.33p one readily 
verifies that p is Tf-invariant. p is regular by observation. The weak continuity of the dynamics in 
the GNS-representation of p will follow from the continuity of the functions of the form 

(3.48) t ^ p{W{gi)W{Ttf)W{g2)), for 51,92,/ G V. 

When Lo > 0, this continuity can be easily observed from the following expresion: 
p{W{gi)W{Ttf)W{g2)) =e^''^s^^92)l^^^^{T^^m~9^)l'i^ 

Note that Tt is differentiable with bounded derivative and that both U and V are bounded. This 
establishes the continuity in the case that a; > 0. 
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As discussed in the introduction of the section, the W*-dynamical system is now defined by 
considering the GNS representation Tip of p. This yields a von Neumann algebra M = 7rp(W(P)). 
The invariance of p implies that the dynamics is implementable by unitaries Ut, i.e., 

(3.50) MniWif))) = U:7rp{Wif))Ut . 

Using Ut, the dynamics can be extended to M.. As a consequence of (|3.48|) . this extended dynamics 
is weakly continuous. 

3.4.4. The case of to = 0. We now discuss the case to = 0. Here, the maps Tt are defined using 
the convolution formula (I3.4ip . By Lemma [3.H Tt is well-defined as a transformation of £^(Z'^), for 
1 < P < 2. Both the group property of Tt and the invariance of the symplectic form a follow in the 
limit w — > by dominated congervence which is justified by Lemma 13. 1[ This demonstrates that 
the dynamics is well defined. 

We represent the dynamics in a state p is defined by (I3.47p . but with the understanding that 
II (f7* — V^*)/|| may take on the value +oo, in which case p{W{f)) = 0. p is still clearly regular. It 
remains to show that the dynamics is weakly continuous. 

Observe that 

(3.51) /• 

follows from (j3.4ip . Using Young's inequality and Corollarv 13.21 it is clear that HTj/ — /|| — > as 
t ^ for any / G l^ill^) with 1 < p < 2. A calculation shows that 

(3.52) {U* - V*){Ttf -f) = Fi* (i/f ) - 5o) - F2 * H^T^^ - iF^ * , 
where 

(3.53) , , 

F2 = .F-^M^^Rel/] , and F3 = M^-,/2J^m[f] . 

A similar argument to what is given above now implies that \\{U* — V*){Ttf — f)\\ — > as t — > 0, 
for any f G T>q, where 

(3.54) Po = {/ G ■■ J^^^M^-^/2J'Re[f] £ £\Z'^)^ . 

No additional assumption on Im[/] is necessary since F^ is convolved with H^^\ Given the form 
of (j3.49p . this suffices to prove weak continuity. In fact, one can check that Tt leaves T>q invariant 
and that if / G Vq, then {U* - V*)Ttf G e'^{Z'^) for all t G R. This estabhshes weak continuity of 
the dynamics, defined on W{T>o). 

Remark 3.3. We observe that, when cj = 0, the finite volume Hamiltonian H'^ ()3.12p is translation 
invariant and commutes with the total momentum operator Pq (see (j3.16p ). In fact, ffj^ can be 
written as 

Hi = PS+ Yl PkPk + l\k)QlQk 

fcGA*\{0} 

fceA2\{o} 

where we used the notation (|3.16p and, for A; / 0, we introduced the operators bk, b\ as in (j3.19p . In 
this case, the operator does not have eigenvectors; its spectrum is purely continuous. By a uni- 
tary transformation, the Hilbert space H-Ki^ (see (I3.13P ) can be mapped into the space L^(]R, dPo! Wfe) 
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of square integrable functions of Pq £ M, with values in Tib- Here, Tib denotes the Fock space gener- 
ated by ah creation and annihilation operators 6^, bk with k ^ 0. It is then easy to construct vectors 
which minimize the energy by a given distribution of the total momentum; for an arbitrary (com- 
plex valued) / G L^{R) with ||/|| = 1, we define ipj £ L^(R,dPo;nb) by setting V/(^'o) = fiPo)^ 
(where Q is the Fock vacuum in Tib)- These vectors are not invariant with respect to the time 
evolution. It is simple to check that the Schrodinger evolution of "0/ is given by e"'^^*^'/ = V'/t 
with ft{Po) = e~**^o^/(Po) is the free evolution of /. In particular, for a; = 0, does not have 
a ground state in the traditional sense of an eigenvector. For this reason, when w = 0, it is not 
a priori clear what the natural choice of state should be. As is discussed above, one possibility is 
to consider first uj ^ and then take the limit a; — > 0. This yields a ground state for the infinite 
system with vanishing center of mass momentum of the oscillators. By considering non-zero values 
for the center of mass momentum, one can also define other states with similar properties. 

3.4.5. Some final comments. The analysis in the following sections and our main result is not 
limited to the class of examples we discussed above. E.g., harmonic systems defined on more 
general graphs, such as the ones considered in (HJ [7] can also be treated. Also note that our choice 
of time-invariant state, while natural, is by no means the only possible. Instead of the vacuum 
state defined in (j3.47p . equilibrium states at positive temperatures could be used in exactly the 
same way. It would also make sense to study the convergence of the equilibrium or ground states 
for the perturbed dynamics and to consider the dynamics in the representation of the limiting 
infinite-system state, but we have not studied this situation and will not discuss it in this paper. 

4. Perturbing the Harmonic Dynamics 

In this section, we will discuss finite volume perturbations of the infinite volume harmonic dynam- 
ics which we defined in Section [3j To begin, we recall a fundamental result about perturbations 
of quantum dynamics defined by adding a bounded term to the generator. This is a version of 
what is usually known as the Dyson or Duhamel expansion. The following statement summarizes 
Proposition 5.4.1 of [1]. 

Proposition 4.1. Let {Ai, at} be a W* -dynamical system and let 5 denote the infinitesimal gener- 
ator of at- Given any P = P* £ Ai, set 5p to be the bounded derivation with domain D{6p) = Ai 
satisfying 6p{A) = i[P,A\ for all A £ M. It follows that 6 -\- 5p generates a one-parameter group 
of *- automorphisms of M. which is the unique solution of the integral equation 

(4.1) af{A)=at{A)+i T af ([P, at_,(^)]) ds . 

Jo 

In addition, the estimate 

(4.2) ||af (A) - at{A)\\ < (el'IH^H - l) P|| 
holds for allt eR and A £ M. 

Since the initial dynamics at is assumed weakly continuous, the norm estimate (14. 2p can be used 
to show that the perturbed dynamics is also weakly continuous. Hence, for each P = P* £ M the 
pair {A4,a^} is also a M^*-dynamical system. Thus, if Pi = P* £ A4 for i = 1,2, then one can 
define af^"*"^^ iteratively. 

4.1. A Lieb- Robinson bound for on-site perturbations. In this section we will consider 
perturbations of the harmonic dynamics defined in Section [3l Recall that our general assumptions 
for the harmonic dynamics on F are as follows. 
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We assume that the harmonic dynamics, t^, is defined on a Weyl algebra W{'D) where V is a 
subspace of ^^(r). In fact, we assume there exists a group Tt of real-hnear transformations which 
leave 2? invariant and satisfy 

(4.3) rfiWif)) = WiTtf) for all /G P. 

In addition, we assume that this harmonic dynamics satisfies a Lieb-Robinson bound. Specifically, 
we suppose that there exists a number ao > for which given any < a < ao, there are positive 
numbers Ca and Va for which 

<cae''»l*l \fix)\\giy)\Fa{dix,y)) 

here the spatial decay in T is described by the function Fa as introduced in Section [2j As we 
discussed in Section [3l the estimate ()4.4p immediately implies the Lieb-Robinson bound 

(4.5) \\[T^{W{f)),W{g)]\\<Cae-^\'\ ^ |/(x)| |5(y)| (d(x, y)) . 

x,yer 

Finally, we assume that we have represented this harmonic dynamics in a regular and r^'^-invariant 
state p for which the pair {Ai,Tf}, with Ai = 7rp{W{'D)), is a I^*-dynamical system. 

Our first estimate involves perturbations defined as finite sums of on-site terms. More specifically, 
the perturbations we consider are defined as follows. 

To each site x S F, we will associate a finite measure px on C, and an element Px G W{'D) which 
has the form 

(4.6) Px= [ W{z6x)fix{dz). 

We require that each fix is even, i.e. invariant under z —z, to ensure self-adjointness, i.e. 
P* = Px- Our Lieb-Robinson bounds hold under the additional assumption that the second moment 
is uniformly bounded, i.e. 

(4.7) sup / \ z\'^ \fix\{dz) < oo . 

xgt Jc 

We use Proposition 14.11 to define the perturbed dynamics. Fix a finite set A C F. Set 

(4.8) P^ = ^P,, 

x6A 

and note that (P^)* = P^ G W('E'). We will denote by t^^^ the dynamics that results from 
applying Proposition 14.11 to the M/^* -dynamical system {A4,Tf} and P^. 
Before we begin the proof of our estimate, we discuss two examples. 

Example. 1) Let fix be supported on [— 7r,7r) and absolutely continuous with respect to Lebesgue 
measure, i.e. fix{dz) = Vx{z)dz. If Vx is in L^([— vr, vr)), then Px is proportional to an operator of 
multiplication by the inverse Fourier transform of Vx- Moreover, since the support of fix is real, Px 
corresponds to multiplication by a function depending only on Qx- 

Example. 2) Let fix have finite support, e.g., take supp(//a;) = {z, —z} for some number z = a + ij3 G 
C. Then 

(4.9) Px = W{z5x) + W{-z5x) = 2 cos(ag^ + ^px) . 
We now state our first result. 



(4.4) 
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Theorem 4.2. Let be a harmonic dynamics defined on T as described above. Suppose that 
(4.10) 



K = sup / \z\'^\fj,x\{dz) < oo , 
xer Jc 



and define the perturbed dynamics r^^^ as indicated above. For every < a < oq, there exist positive 
numbers Ca and Va for which the estimate 



(4.11) 



holds for all t gM. and for any functions f,g £ T>. 



Here the numbers Ca and Va are as in (j4.4p . whereas Ca is the convolution constant as defined in 
2.2p with respect to the function Fa- 



Proof. Fix t > and define the function : [0,t] W(P) by setting 



(4.12) 



"^tis) 



(A) f^o 



t-s 



{W{f))),W{g) 



It is clear that interpolates between the commutator associated with the original harmonic 
dynamics, at s = 0, and that of the perturbed dynamics, rj^^^ at s = t. A calculation shows that 



(4.13) 



5][rf)([P.,VF(W)]),VF(g) 



where differentiability is guaranteed by the results of Proposition 14.11 The inner commutator can 
be expressed as 



[P^,W{Tt^sf)] = [ [W{z6,),W{Tt^sf)]l^x{dz 
Jc 

= W{Tt.sf)Ct-sM)- 



(4.14) 
where 

(4.15) CU-M) = ^t-s-M) = ^ W[z5^) |e-m-./,.<5.) _ i| ^^(^^) e ■^^(p) 

Thus satisfies 



(4.16) 



d_ 
ds 



^,(s)=zj^^t(.)ri^) {Ct-s-Af)) 



The first term above is norm preserving. In fact, define a unitary evolution Ut{-) by setting 
"^-Utis) = -iYl ^i""^ i^t-sAf)) Ut{s) with Ut{0) = 11 . 



(4.17) 



ds 



It is easy to see that 

(4.18) 4:. iMs)Ut{s)) = . ^i^^ (WiTt.sf)) [rf ) {CtsAf)) , ^^(5)1 Ut{s) 



ds 



x£A 



and therefore. 



(4.19) ^if,{t)Ut{t) = MO) + iY ) (W^(ri_,/)) )(£*_,;,(/)), VF(<7) Utis) ds . 

x6A-^o ^ ^ 
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Estimating in norm, we find that 



.(A) 



(4.20) 



iW{f)),W{g)]\\ <\\[THW{f)),W{g) 



+ E 

x6A 



rf) (£*_,;.(/)), VF(g) 



ds . 



Moreover, using (j4.15p and the bound (j4.4p . it is clear that 



ri^) {Ct-sAf)),W{g) 



Va{t-s) 



Y,\fix')\Fa{d{x,x')) 



(4.21) 



x'er 



[rf) {W{z6^)),W{g) 
holds. Combining (j4.2ip . (j4.20p . and (j4.5p . we have proven that 

rf ) {W{f)) , W{g)] < c,e^"* Yl 1^(^)1 \9iy)\ Fa {d{x, y)) 



\fj,x\{dz) 



x,y 



(4.22) 



+ Ca J^|/(x')l5]Fa(d(rE,rE')) / e-"( 
x'er xeA 

[rf) (VF(^4)),iy(g)] |||//,|(a!z)ds. 



Following the iteration scheme applied in [TO], one arrives at (|4.1ip as claimed. 



□ 



4.2. Multiple Site Anharmonicities. In this section, we will prove that Lieb-Robinson bounds, 
similar to those in Theorem 14.21 also hold for perturbations involving short range interations. We 
introduce these as follows. 

For each finite subset X C F, we associate a finite measure i^x on C"''- and an element Px £ W{V) 
with the form 



(4.23) Px= / Wiz-6x)f^x{dz), 
where, for each z € C"'*', the function z • (5x : F — > C is given by 

(4.24) {z-6x)ix)=Yl 



x'&X 



Zx if X G X, 
otherwise. 



We will again require that ^x is invariant with respect to z 
In analogy to (14. Sp . for any finite subset A C F, we will set 



-z, and hence, Px is self- adjoint. 



(4.25) 



P' 



XcA 



where the sum is over all subsets of A. Here we will again let r^^^^ denote the dynamics resulting 
from Proposition 14 . 1 1 applied to the VF*-dynamical system {}A,t^^ and the perturbation P^ defined 
by (105]) . 

The main assumption on these multi-site perturbations follows. There exists a number ai > 
such that for all < a < ai, there is a number > for which given any pair xi,X2 G F, 



(4.26) 



E 

XCT: 
Xl,X2&X 



\Zxi\\Zx2\\l^x\{dz) < KaFa {d{xi,X2)) 
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Theorem 4.3. Let be a harmonic dynamics defined on T . Assmue that |^.g6] j holds, and that 
T^^^ denotes the corresponding perturbed dynamics. For every < a < min(ao,ai), there exist 
positive numbers Ca and Va for which the estimate 



(4.27) 



x,y 



holds for all t and for any functions f,g £ V. 

The proof of this result closely follows that of Theorem 14.21 and so we only comment on the 
differences. 



Proof. For f,g £ T> and t > 0, define ■ [0, — > W{T>) as in (|4.12|) . The derivative calculation 
beginning with ()4.13p proceeds as before. Here 



Ct-sMf) = I W{z- 5x) \ e-(^*-/'-^^) - 1 



is also self-adjoint. The norm estimate 



.(A) 



(4.29) 



{W{f)),W{g)\\\ <\\[Tj^{W{f)),W{g) 



+ E / \r^s''H^t-sMf)),Wig) 



ds . 



holds similarly. With (|4.28|) . it is easy to see that the integrand in (|4.29|) is bounded by 



(4.30) c,e'"'(*-^) Yl 1/(^)1 E {d{x,x')) ! 



T'^''Hw{z-5x)),W{g) \fix\{dz), 



the analogue of (|4.21|) . for < a < qq. Moreover, if < a < min(ao,ai), then 



(4.31) 



rf^(W(/)),W^(9)]||<Cae^"'* E \f{x)\\9{y)\Fa{d{x,y)) 

+ 1/(^)1 E E^4rf(^'^'))x 

x&r xcAx'(^x 

X Te-^t-) f \4''){W{z.6x)),W{g) 
Jo Jcx " L 



\fix\{dz) ds . 
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The estimate claimed in (|4.27p follows by iteration. In fact, the first term in the iteration is bounded 
by 

(4.32) 

ft 



X / kxil (cae""" 'Y\z^^\\g{y)\Fa{d{x2,y))^\iix\{dz)ds 

■^'^^ X2GA- y 

<Cat-Ca^'-'Y}f^^)\\9{y)\ E Fa{d{x,Xi))Fa{d{x2,y)) E j Jz^,\\Zx,Ux\{dz) 

x,y xi,x2er xcr-. 

x\ ,X2£X 

<^€aCat-Cae^"*^|/(a;)||5(y)| J2 Fa{d{x,xi))Fa{d{xuX2))Fa{d{x2,y)) 
x,y xi,x2eT 



< 

x,y 



The higher order iterates are treated similarly. □ 



5. Existence of the Dynamics 

In this section, we demonstrate that the finite volume dynamics analyzed in the previous section 
converge to a limiting dynamics as the volume A on which the perturbation is defined tends to T. 
We state this as Theorem 15.11 below. 



Theorem 5.1. Let be a harmonic dynamics defined on W{i^{T)) as described in Section \4- 1\ 
Let {An} denote a non- decreasing, exhaustive sequence of finite subsets ofT. Consider a family of 
perturbations P^" as defined in \4.25 ) and 1^4.23 ) which satisfy \4.26 ). Suppose in addition that 



(5.1) M = SUpV/ \Z:,\\llx\{dz) < CO . 

xer Jc^ 



x&X 



Then, for each f G ^^{T) and t G M fixed, the limit 

(5.2) hm {W{f)) 

n — >oo 

exists in norm. The limiting dynamics, which we denote by Tt, is weakly continuous. 

It is important to note that since the estimates in Theorem 14. 31 are independent of A, the limiting 
dynamics also satisfies a Lieb- Robinson bound as in (14.270 . We now prove Theorem 15.11 

Proof. Fix a Weyl operator W{f) with / G ^^(T). Let T > and take m < n. Iteratively applying 
Proposition 14.11 we have that 



(5.3) 



^f")(Ty(/)) = rf-)(Ty(/)) + z^Vi^") ([p^"\^-,Ti!^-\wif)) 



ds . 
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for all —T < t < T. The bound 



(5.4) 



< 



E 



XCA„: 

XnA„\A™^0 



< Cae 



(Va+CaKaC^)(t—s) 



< Cae 



iVa+CaKaC^)(t-s) 



Wiz-5x),ri!^T\wif))\ \fix\idz) 
XnA„\A„7^0 

El/(^)l E i^a((i(x,y))5] / 



m 



xcr: 



xer j/eA„\Am 

follows readily from Theorem 14.31 and assumption (jS.ip . For / G ^^(T) and fixed i, the upper 
estimate above goes to zero as n, m — > oo. In fact, the convergence is uniform for t G [— T, T]. This 
proves (|5.2p . 

By an e/3 argument, similar to what is done at the end of Section [21 weak continuity follows since 
we know it holds for the finite volume dynamics. This completes the proof of Theorem 15. 1[ □ 
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